Phase diagram for a Bose-Einstein condensate moving in an optical lattice 



O 
O 

<N 

a 



5-H 



as 



a 
o 
o 



> 
m 

O 
l> 
O 



X 
S3 



Jongchul Mun, Patrick Medley, Gretchen K. CampbellQ Luis G. MarcassaQ David E. Pritchard, and Wolfgang Ketterle 

MIT-Harvard Center for Ultracold Atoms, Research Laboratory of Electronics, 
and Department of Physics, MIT, Cambridge, Massachusetts 02139, USA. 
(Dated: February 1, 2008) 

The stability of superfluid currents in a system of ultracold bosons was studied using a moving optical lattice. 
Superfluid currents in a very weak lattice become unstable when their momentum exceeds 0.5 recoil momentum. 
Superfluidity vanishes already for zero momentum as the lattice deep reaches the Mott insulator(MI) phase 
transition. We study the phase diagram for the disappearance of superfluidity as a function of momentum and 
lattice depth between these two limits. Our phase boundary extrapolates to the critical lattice depth for the 
superfluid-to-MI transition with 2% precision. When a one-dimensional gas was loaded into a moving optical 
lattice a sudden broadening of the transition between stable and unstable phases was observed. 
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The realization of condensed matter systems using ultra- 
cold atoms brings the precision and control of atomic physics 
to the study of many-body physics. Many studies have fo- 
cussed on Mott insulator physics, an important paradigm for 
the suppression of transport by particle correlations. Previous 
studies of the superfluid(SF)-to-Mott insulator(MI) transition 
in optical lattices with ultracold bosons HI 0, 
adressed the quenching of superfluidity below a critical lattice 
depth. Here we extend these studies into a second dimension 
by studying stability of superfluid current as a function of mo- 
mentum and lattice depth as suggested in ref. 190 . These trans- 
port measurements show the stability of the quantum phase in 
a moving system, which is far from equilibrium. 

Transport measurements extend previous work on station- 
ary systems in two regards. First, superfluidity near the MI 
transition has only been indirectly inferred from coherence 
measurements, whereas in this work, we characterize the su- 
perfluid regime by observing a critical current for superfluid 
flow. Second, previous studies lll|2l|3i|^l^lal2llSl] were not 
able to precisely locate the phase transition, since the observed 
excitation spectrum and atomic interference pattern did not 
abruptly change |3, 5, 6], partially due to the inhomogeneous 
density. In contrast, the sudden onset of dissipation provides 
a clear distinction between the two quantum phases. In the 
SF phase, current flows without dissipation if the momentum 
does not exceed a critical momentum, while in the MI phase 
the critical momentum vanishes and transport is dissipative. 

Bosonic atoms in an optical lattice are often described by 
the Bose-Hubbard Model where the tunneling between near- 
est neighbour lattice sites is characterized by the hopping ma- 
trix element J and the repulsive interactions by the on-site 
matrix element U JH QH Qjl [12] • The dimensionless inter- 
action energy u = U/ J determines the quantum phase of the 
system. For u > u c , the system is in the MI phase; for u < u c , 
it is in the SF phase. 

For weak interactions (u —> 0), the system approaches 
single-particle physics in a periodic potential well described 
by Bloch states and band structure. The critical momentum 
for a stable current-carrying state is 0.5 p r (p r = h/X is the 
recoil momentum of an atom, where A is the wavelength of 
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FIG. 1 : Phase diagram showing the stability of superfluid flow in an 
optical lattice and the experimental procedure. The grey curve shows 
the predicted boundary between superfluid flow and dissipative flow 
phases for a three-dimensional gas with a commensurate filling of 
N = 1 atom per site | 9]. The solid (dashed) arrows illustrate the 
experimental trajectory used for small (large) lattice depths (see text 
for details). 



the optical lattice light) 11311 . At the critical momentum, it be- 
comes possible for two atoms in the same initial Bloch state 
to scatter into two other states and conserve energy and quasi- 
momentum Il4i fish . For sufficiently deep lattices this oc- 
curs when the effective mass becomes negative. Instabilities 
in a ID optical lattice were studied theoretically using a linear 
stability analysis of the Gross-Pitaevskii Equation EE] , 
and experimentally rfvL fTyh . The theoretical studies predicted 
that for increasing lattice depth or increasing atomic interac- 
tions the stability of superfluid flow should increase 
the dynamic instability would stay near 0.5 p r whereas the 
Landau critical velocity and therefore the energetic instability 
would shift to larger momenta (For more discussions on dy- 
namic and energetic instability, see ref. Qji]). However, 
these analyses neglect the growing importance of quantum 
correlations for larger lattice depth which lead to the SF-MI 
phase transition, where the critical momentum for a superfluid 
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current vanishes. In this paper, we study the decrease of the 
critical momentum from its value for the weakly interacting 
regime towards zero at the MI transition. (Fig. [T]) 

Most studies of the SF-MI phase transition monitor the co- 
herence in the superfluid phase through an interference pattern 
observed in the ballistic expansion resulting from a sudden 
turn-off of the confining potential and lattice. Previous obser- 
vations of the phase transition found the experimental transi- 
tion point to lie in the range between 10 and 13 Er (with the 
recoil energy defined as Er = p%/ 2m where m is the atomic 
mass) 0]. This uncertainty is related to the inhomogeneous 
density profile of trapped atoms and to the fact that the visi- 
bility of the interference extends beyond the transition point 
due to short-range coherence in the MI phase |6] . It has been 
suggested that observed kinks in the visibility are linked to the 
formation of the MI shells with occupation numbers N = 2 
and 3 1 6]. Several authors have suggested other features in the 
momentum distribution beyond coherent interference peaks as 
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a more distinct signature of the phase transition 12QLI21H . Here 
we show that the disappearance of the critical momentum for 
superfluid flow provides such a signature and allows the de- 
termination of the transition point with high precision. 

Our measurement was not limited by the inhomogeneous 
density profile. For our range of lattice depths, low critical 
momenta and the onset of dissipation occur only near the for- 
mation of MI shells with integer occupation numbers N |9|]. 
The onset of dissipation related to the N = 1 domains occurs 
at smaller momentum than for other N domains. For instance, 
with increasing momentum p the N = 1 domain becomes un- 
stable first, and this triggers dissipation over the whole atomic 
cloud [9] . Therefore, the breakdown of superfluid flow in the 
system was determined by the formation of the N = 1 do- 
main and was not smeared out by the inhomogeneous density. 
Our criterion, the sudden onset of dissipation, depended on 
the formation of an insulating shell surrounded by a superfluid 
region, which occurs only in the inhomogeneous case. 

In our experimental setup, a Bose-Einstein condensate of 
87 Rb atoms in the 5Si/ 2 1) state was prepared and 
trapped in a combination of an Ioffe-Pritchard magnetic trap 
and an optical dipole trap. The number of atoms in the 
BEC was typically 2 x 10 5 , resulting in a maximum filling 
factor N of around 3. The magnetic trap frequencies were 
u x ,y = 40 Hz radially and oj z = 4.6 Hz axially. The laser 
beam for the optical dipole trap was oriented along the x- 
axis. This laser beam was retroreflected and the polariza- 
tion of the retroreflected beam was rotated in order to mini- 
mize interference between the two beams. Along the verti- 
cal direction (y-axis) a lattice was formed by a retroreflected 
laser beam. For the z-axis, a moving lattice was created by 
introducing a small frequency detuning Sf between the two 
counter-propagating laser beams using acousto-optical modu- 
lators driven by phase-locked frequency generators. The 3D 
optical lattice was ramped up in the following way: For the 
lattice along the x-axis, the polarization of the retroreflected 
beam was rotated to increase the interference contrast. For the 
other two axes, the power was increased exponentially in 160 
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FIG. 2: Determination of the critical momentum of superfluid flow. 
Shown is the condensate fraction as a function of a momentum p. 
(a) Condensate fraction with u/u c — 0.61 for a variable number 
of cycles of the momentum modulation (one cycle: cross marks and 
blue line, two cycles: squares and purple line, three cycles: triangles 
and red line). A dashed vertical line indicates the critical momen- 
tum where instability begins to occur. The two and three-cycle data 
are offset vertically for clarity. These data were fitted with an error 
function to guide the eye. (b) Images of interference patterns released 
from an optical lattice at u/u c = 0.61 moving with variable momen- 
tum. Instability occurred between p = 0.31p r and 0.32p r . Some of 
the triangular data points in (a) were obtained from these images, (c) 
Condensate fraction on a log-log scale for two different interaction 
strengths. 



ms. All lattice beams were derived from the same laser oper- 
ating at A = 1064 nm and had an 1/e 2 waist of 100 ~ 200 
pm. The lattice depth was calibrated with 1 % accuracy by 
applying a 12.5 ps lattice laser pulse to a BEC and compar- 
ing the observed Kapitza-Dirac diffraction pattern of a BEC 
to theory. Fits of more than twenty diffraction patterns for the 
lattice depth from < 1 Er to 40 Er resulted in an accuracy of 
the calibration of 1 %. 

For transport measurements, we moved an optical lattice 
1 17l I22H which provides more flexibility to change the mo- 
mentum than exciting a dipole oscillation by displacement of 
the BEC A moving optical lattice with velocity v = 

A • 5f/2 was created along the long axis of the BEC by intro- 
ducing a small frequency detuning Sf between two counter- 
propagating lattice beams. If the velocity v(t) changes slowly 
enough not to induce interband excitations, the initial Bloch 
state \p = 0) of the condensate in the optical lattice adiabati- 
cally evolves into the current-carrying state \p(t) = —mv(t)) 
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where p is the quasi-momentum. For increasing lattice depth, 
the effective mass of atoms m* = [d 2 E(p)/dp 2 ]~ 1 increases, 
and the group velocity v g = —(m/m*)v(t) decreases. As a 
consequence, atoms prepared in a moving lattice with quasi- 
momentum p — —mv travel in the frame of the moving lattice 
with v g and in the lab frame with velocity Av = v ± v g = 
(1 — m/m*)v, which approaches v in a deep lattice. Con- 
sequently, we observed that in a deep moving lattice atoms 
were dragged along to the edge of the trapping region limiting 
the experimental time scale to probe for dissipation. This be- 
came an issue for larger values of p and was addressed by first 
ramping up the lattice with p = and then alternating the ve- 
locity of the moving lattice, thus performing a low-frequency 
AC transport measurement instead of DC. 

We have used two sets of experimental procedures (Fig. [T]), 
and our results were consistent for both. Close to the SF-MI 
phase transition, the lattice was increased to Vi aU with a fixed 
(and small) value of momentum p (dahsed arrows in Fig. Q]). 
After a variable hold time thoid at Viatt the lattice was ramped 
down to zero, and the magnetic trap switched off. After 33 ms 
of ballistic expansion, the atoms were imaged and the conden- 
sate fraction was determined as a function of momentum. For 
smaller lattice depths, the lattice was ramped up with p = 
(Fig ID- Then a sinusoidal momentum modulation of the mov- 
ing lattice with amplitude pm was applied by modulating the 
frequency detuning Sf between the counter-propagating lat- 
tice beams. The 10 ms period of this momentum modulation 
was slow enough to meet the adiabaticity condition, but fast 
enough to limit the displacement of the atomic cloud to less 
than a few fim. Both the trapping potential and the optical 
lattice were then turned off suddenly. After 33 ms of ballis- 
tic expansion, the condensate fraction of the center peak of 
the superfluid interference pattern was recorded as a function 
of the momentum modulation amplitude pw Several cycles 
(typically, three to five) of the momentum modulation were 
applied to obtain a high contrast between the stable and dissi- 
pative regimes (Fig. [2(a)). 

Fig. [2(a) shows how the transition between superfluid and 
dissipative currents became sharper with increasing number of 
cycles of the momentum modulation. The critical momentum 
was determined from a log-log plot of the condensate frac- 
tion as a function of momentum p (Fig. [21 (c)). For low p, 
the condensate fraction is constant. For high p, the conden- 
sate fraction decreases. The intersection between two linear 
fit functions was taken as the critical momentum. Our result 
was found to be independent of the time period and number 
of cycles of the momentum modulation at a few percent level. 

In the MI phase, stable superfluid flow is not possible and 
the critical momentum should vanish. However, using the 
procedure described above, we measured a small critical mo- 
mentum of 0.02 p r for lattice depths Vi a tt = 14, 15, 16 Er. 
Up to this momentum, the SF-MI phase transition remained 
reversible. We interpret the non-zero critical momentum as 
a finite-size effect. For our cloud size of 60 fim, the cor- 
responding Heisenberg momentum uncertainty of 0.018 p r 
agrees with our measured critical momentum. In cold atom 
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FIG. 3: Critical momentum for a condensate in a 3D lattice. The 
solid line shows the theoretical prediction for the superfluid region. 
The horizontal solid line is a fit to the data points in the MI phase. 
(Inset) Fit of critical momenta near the SF-MI phase transition. 



experiments, some sloshing motion of the cloud in the trap- 
ping potential is unavoidable. The momentum uncertainty de- 
termined above indicates how much sloshing motion can be 
tolerated without affecting the observed phase transition. 

The critical lattice depth for the SF-MI phase transition 
can be determined as the point where the critical momentum 
vanishes. Using the predicted functional form [9] of the ap- 
proach towards zero, p c ex yl — u/u c , as a fit function for 
the data points close to the SF-MI phase transition (the data 
points shown in the inset of Fig. [5]) we determined the critical 
value u c = 34.2 (±2.0) corresponding to a lattice depth of 
13. 5 (±0.2) Er B29ll . Our result agrees with the mean field 
theory prediction u c = 5.8 x 6 = 34.8 for N = 1 SF- 
MI phase transition [1] and deviates by 2 a from the pre- 
dictions of Uc = 29.34(2) of quantum Monte Carlo(QMC) 
simulation |26|], which includes corrections beyond the 
mean field theory. This demonstrates that our method has the 
precision to identify non-mean field corrections. However, to 
turn precision into accuracy, experiments or QMC simulations 
[HclEUE] have to address corrections due to finite size, finite 
temperature, and finite time to probe the onset of the instabil- 
ity ll26n . In our experiment, these corrections seemed to be 
small, but have not been characterized at the level of 1% in 
lattice depth. 

The mean-field prediction for stable superfluid flow in ID 
is similar to that for the 3D system |9[]. However, it is well 
known that fluctuations play a much more important role in 
ID. For studying a ID system, we prepared an array of one- 
dimensional gas tubes by ramping two pairs of optical lattice 
beams up to lattice depths of V x = V y = 30 Er suppressing 
hopping between the tubes. After a hold time of 10 ms, a mov- 
ing optical lattice was ramped up along the z-axis. As in our 
3D experiment, a momentum modulation was applied, after 
which the moving optical lattice was ramped down to zero, 
followed by the other two optical lattices. The condensate 
fraction was determined after 33 ms of ballistic expansion as 
a function of the momentum modulation amplitude. The crit- 
ical momentum, where the onset of dissipation begins, was 
identified from a log-log plot as in the 3D case. Since the 
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FIG. 4: Critical momentum for a ID gas in an optical lattice, (a) 
The grey line indicates the mean-field theory prediction. The in- 
teraction strengths are normalized by the mean-field prediction for 
u c = 5.8 x 2 |liyt|. Squares(crosses) represent the measured crit- 
ical momentum (the center of the transition). Measurements were 
taken at lattice depths of 0.25, 0.50, 0.75, 1.0, 2.0 Er. The lines 
between crosses and squares indicate the width of the transition re- 
gion, (b) Condensate fraction measured at 0.25 Er (left) and 0.75 
Er (right). The data were fitted with an error function. Squares (the 
critical momentum) and crosses (the center of the transition) are in- 
dicated on the plots. 

transitions became very broad, we characterized them by an 
error function fit, with the center of the fitted error function 
taken as the center of the transition (Fig. 0]). 

In the ID system, at a very shallow lattice depth of 0.25 
Er (corresponding to u/u c = 0.08) a sharp transition was 
observed, and the measured critical momentum agreed very 
well with the prediction |9l [27|] of a critical momentum of 
0.39 p r . However, a slight increase of the interaction strength 
(to u/u c = 0.09 at a lattice depth of 0.5 Er) led to a signifi- 
cant decrease of the critical momentum as well as a dramatic 
broadening of the transition as shown in Fig. |4l For lattice 
depths larger than 2 Er, the transition became very broad and 
showed complex behavior, and we could not obtain quanti- 
tative fits. Our results show a significant deviation from the 
mean-field theory predictions and are in agreement with pre- 
vious experiments 12411 where damped dipole oscillations of a 
ID Bose gas in an optical lattice were observed, and the damp- 
ing grew rapidly even at very shallow lattice depths around 
0.25 E R . 

The observed broadening of the transition confirms theo- 
retical studies which emphasize the importance of quantum 
fluctuations in the ID system. Quantum tunneling out of 
metastable states which are ignored in the mean-field descrip- 
tion can lead to a decay of the superfluid current at very low 
momenta [27]. In addition to quantum fluctuations, thermal 
fluctuations provide a mechanism for current decay |27|]. In 



our experiment, we used a "pure" BEC without a discernible 
thermal component. The close agreement with T = predic- 
tions indicates that thermal fluctuations were not dominant. 

In conclusion we have used transport studies to connect 
a well-known dynamical instability for weakly interacting 
bosons with the equilibrium superfluid to Mott insulator tran- 
sition. A comparison of 3D and ID systems confirms the ap- 
plicability of a mean-field description in 3D and the crucial 
importance of fluctuations in ID. The disappearance of super- 
fluid currents at the SF-MI phase transition precisely located 
the phase transition. Our results illustrate the control and pre- 
cision of condensed matter physics experiments done with ul- 
tracold atoms and their suitability to test many-body theories. 
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